o 
o 

(N 
o 

Q 



in 



> 

o 
o 

o 



The Time Periodic Solution of the Burgers Equation 
on the Half- Line and an Application to Steady Streaming 

A.S. Fokas^ and J.T. Stuart^ 
^Department of Applied Mathematics and Theoretical Physics 
University of Cambridge, 
Cambridge, CB30WA, UK 
t.fokas@damtp.cam.ac.uk 
^Department of Mathematics 
Imperial College 
London, SW72BZ, UK 
t . st uar t @ i mp erial.ac.uk 

October 2004 

Abstract 

The phenomenon of steady streaming, or acoustic streaming, is an important physical phenomenon studied 
extensively in the literature. Its mathematical formulation involves the Navier-Stokes equations, thus due to the 
complexity of these equations is usually studied heuristically using formal perturbation expansions. It turns out 
that the Burgers equation formulated on the half line provides a simple model of the above phenomenon. The 
physical situation corresponds to the solution of the Dirichlet problem on the half-line, which decays as a; ^ oo and 
which is time •periodic. We show that the Dirichlet problem, where the usual prescription of the initial condition is 
now replaced by the requirement of the time periodicity, yields a well posed problem. Furthermore, we show that 
the solution of this problem tends to the "inner" and "outer" solutions obtained by the perturbation expansions. 



1 Introduction 

\^ • The phenomenon of "steady streaming", or "acoustic streaming" as it is sometimes called, dates back to the 19th 
. ^ ^ century with the work of Kundt on the circulations of air in tubes (see Rayleigh 1883). There are related phenomena 
of circulations in the flows under water waves with a free surface (Longuet-Higgins 1953, 1960; Hunt and Johns 

^ \ 1963). Following Rayleigh, Schlichting (1932) did both experimental and theoretical work on a particular case of the 

- - ' flow induced by a circular cylinder oscillating along a diameter. 

The nature of the phenomenon, which is common to all the papers mentioned above, may be explained as follows. 
If the flow of fluid (liquid or gas) is periodic in time but with zero mean, and is parallel to a rigid surface, then the 
governing equations ensure that a thin periodic boundary layer is formed in the neighborhood of that surface. If the 
time frequency is cj-radians per second and v is the kinematic viscosity then the thickness of this layer is of order 
(ylw)^!"^ . This boundary layer is known as a Stokes layer. A case of particular simplicity is that in which the flow 
is parallel to a flat rigid surface, whose velocity is Uw cos wi parallel to itself, t being the time and IJ^ a constant. 
Then the velocity in the fluid neighboring the surface is 

u = J7iue~'' cos(wt — ?7), (1.1) 

where 

r; = 2;(w/2z.)i/^ (1.2) 

y being the coordinate normal to the surface. This is the prototype Stokes layer (Stokes 1851), see also (Stuart 1963, 
Benney 1964). If a Galilean tranfsormation is imposed such that the surface is at rest but the fluid has velocity 



1 



U COS wt parallel to that surface as y and rj tend to infinity, then the velocity in the fluid is 



U[coswt—e ^ cos{'wt — T])]. (1.3) 

In the above description, the fluid motion is intrinsically linear so that the simple form of the Stokes layer emerges 
exactly: the nonlinear terms are identically zero for (1.1) and (1.3). 

The situation is quite different, however, if the velocity (as y and rj tend to infinity) has the form 

U{^) cos wt, (1.4) 

where ^ is the coordinate parallel to the surface: the latter may be curved and its curvature is neglected if the 

Stokes layer thickness (z//w)^/^ is small compared with a typical length d: thus v/wiP is supposed to be small. The 
nonlinear term in the Navier-Stokes equation, namely u-Vu, where u is the vector velocity, is not zero however. It has 
terms of two types: (i) there are terms proportional to exp(2iwf) and to exp(— 2iwf) which generate corresponding 
flow components; the equations that govern these components do have solutions which tend to zero at the edge of the 
Stokes layer as is required, (ii) a mean term is also generated, which drives the steady streaming; we can refer to this 
mean term, which is an effect of rectification, as a Reynolds stress (or, rather a derivative of the Reynolds stress). 
An equation can be obtained for the mean flow generated by this Reynolds stress but, somewhat paradoxically, it 
is not possible to obtain a solution whose component of velocity that is parallel to the surface, tends to zero as the 
edge of the Stokes layer is approached. Rather the best that can be achieved is to ensure that the solution is finite 
for that component of velocity. 

Indeed it is found that, as the edge of the Stokes layer is approached, the velocity component approaches 

Since ^ has a typical length d, this velocity component has scale UQ/wd, where Uq is the scale of U{£,). Now in 
problems of the Navier-Stokes equations, the famous Reynolds number plays a significant role; it is the product of a 
velocity and a distance, divided by the kinematic viscosity. In the present case we have 

= {U^/wd){d/v) = U^/wv. (1.6) 

This parameter is known as the steady-streaming Reynolds number. Stuart (1963, 1966) showed the importance of 
this concept for the calculation of the flow forced by (1.5) outside the Stokes layer. 

In the papers of Rayleigh (1883) and Schlichting (1932), it is implicit that the parameter Rg is small, and the 
calculation of the flow outside the Stokes layer is performed on that basis. Even so, Schlichting's paper makes it 
clear that he was aware that his theory for small values of Rg is inappropriate as an explanation of his experimental 
work, for which Rg = 250. 

Stuart (1963, 1966) showed how an "outer-boundary-layer theory", which is valid for large values of Rg, can 
be used to obtain a solution for which (1.5) applies at the edge of the Stokes layer but which tends to zero as the 
distance form the surface tends to infinity. Implicit in that work is the idea that this "outer" (or "steady-streaming') 
boundary layer is much thicker than the Stokes layer, indeed by a factor i?y^. Following Stuart's work, Riley (1998) 
has pursued problems of this type in much detail. 

A problem can be posed that shows many of the characteristics that are outlined above, but without the compli- 
cations of the full Navier-Stokes equations. The relevant equation is a Burgers type equation for u{x,t), namely 



ut + (3{u - kP)ua; = -Wxx, (1-7) 

neter which plays 

boundary conditions are 



where /3 is a small parameter which plays a role analogous to Rs above, Rg being large and k is 0(1). The 



X = : u(0, t) = cosL 

a; — > oo : u(oo,f)— +0, (1.8) 



and 

u(x, t) = u{x, t + 27r). 
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The Cole-Hopf transformation 
gives 

with the boundary conditions 



u- kP = -e^/l30 (1.9) 
0t - \o., = 0, (1.10) 



x = Q: e^{0,t) + l3{cost - kl3)e{0,t) = 0, 

x^oo: e^{oo, t)/e{oo, t) = kf3'^, (1.11) 

e^{x,t)/6{x,t) = e^{x,t + 2Tr)/e{x,t + 2Tr). 

The equation (1.10) is simply stated, but the boundary and periodicity conditions are complicated. We return to a 
discussion of this problem in later sections. 

In the meantime we note that a formal solution to (1.7), (1.8) can be obtained by expanding u in a power series 
in 13 

u = uo{x,t) + Pui{x,t) + P'^U2{x,t) H , 

and it is quickly found that 

uq = cos{t - x). (1-12) 
At order /?, the ui terms involving e^'* and e~^** give no difficulty, but the steady (or rectified) term has to satisfy 

the solution which is bounded as a; — > oo and which is zero at a; = is 

u, = \{l-e-'n- (1-13) 

We note that as x ^ oo, u ^ j(3 plus higher order terms. A re-scaUng of (1.7) for the steady part of the solution 
with 

u{x) = I3v{x), X = 

yields 

{v-k)vz=Vzz, (1.14) 

with the boundary conditions 

^^0: v^l/A, 

z^oo: v^O. (1.15) 

The solution of (1.14) subject to (1.15) is 

V = 2k/[{8k-l)e^ + 1]. (1.16) 

We note that the length scale of this "outer" region is (2/3^)^^ times the scale of the "inner" region of (1-12) and 
(1.13). Also z — !■ in (1.15) corresponds to a; ^ oo in (1.13), in the sense of /? ^ with x fixed and f3 ^ with z 
fixed giving an equivalence. 

This heuristic argument will be justified in later sections by solving explicitly Burgers equation (1.7) with the 
conditions (1.8). We note that the mathematical novelty of this problem is that it is posed on the half-line and that 
it requires periodicity in t. We recall that the usual Dirichlet problem for the Burgers equation on the half-line was 
analysed by Calogcro and Dc Lillo. In their formulation, in addition to the Dirichlet boundary condition, one also 
specifies the initial condition M(a;,0). In the present situation, instead of specifying u(x,0), one requires periodicity. 

In section 3 we will analyse the time periodic solution of the Biirgers equation on the half-line with two different 
boundary conditions at the origin: we specify either (a) the integral of u{x,t), or (b) u(0,t). The physical problem 
corresponds to the case (b), however we have also included case (a), because for this case the coefficients An of the 
associated Fourier series (given by the representation (3.3)) can be computed explicitly. For case (b) the coefficients 
A„ satisfy a second order difference equation, see equation (3.4b). We will show that this equation has a unique 
solution by imposing the condition that ^ as n ^ oo (which is needed for the convergence of the series). 
Although we cannot give the explicit form of An in general, we will compute An in the case that /? is small. In this 
case we will show that the associated series converges and that the solution u{x, t) tends to the "inner" and "outer" 
solutions obtained by the perturbation expansion, see equations (5.15) and (5.16). 
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2 From Burgers to the Heat Equation 



Proposition 2.1. Let (3 and c be positive constants. Let the real-valued function u{x,t) solve 

Ut = -jUxx - PuUx + CUx, < X <oo, t>0, 

u ^ as X ^ 00, 
u{x, t) is 27r periodic in t, 
and one of the following boundary conditions: 



or 



/■oo 

/ u{x,t)dx = cost, t>Q, 
Jo 



u{0,t) = cost, t>0. 



Define ip{x, t) by 



^{x,t) = e-^r^^^^^t)di-l. 
Then (p{x, t) is a real-valued function and solves 

= ^<fxx + ct^x, < a; < oo, t>Q, 
— > as a; — > oo, 

if{x, t) is 2-K periodic in t, 
and satisfies one of the following boundary conditions, respectively: 

(/?(0,i) =e'^'=°"*-l, i>0. 



or 



iPx{Q,t) + f3costip{Q,t) +j3cost = 0, t>Q. 



Proof The fiuiction Lp is well defined and <^ — > as a; — > oo. 
If E is defined by 



then 



thus 
Also 

Equation (2.2a) implies 
Furthermore, 



£; = exp -/3 

iPx = -puE, 

Vxx = -PuxE + 0^u^E, 

J oo 

1 

Vt - -j^^xx - C(p, 



J OO 



(^j utdi - ^ux + \(iu^ - ci)j Efi = 0. 



<Px(0, t) = -0u{O, t) exp 



Jo 



(2.1a) 

(2.1b) 
(2.1c) 

(2.2a) 

(2.26) 
(2.3) 

(2.4a) 

(2.4b) 
(2.4c) 

(2.5a) 
(2.56) 
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(3 / uiiM 



hence 



which imphes (2.5b). 



<^(0, t) = exp 



3 The Solution for u{x,t) 

Proposition 3.1. Assume that the positive constants /3 and c, are such that A„ defined by 

A„ = -c+ (c^ + 2in)5, n=l,2,--- 
has 1 value with ReAn < 0. Then the solution of (2.1)-(2.2) is given by 

1 ipx{x,t) 



u{x, t) = — 



I31 + <f{x,t)' 



where 



A „\-nX+int _i_ A p- 



■iX—int^ 



and Ao, {Aj}Q° are defined as follows: 

in the case of the boundary condition (2.2a) 



Ao = -2c; 



JO 



Ar,.= I e"'"* [e'^^"^*- 1] dt, 

whereas in the case of the boundary condition (2.2b) 

2A„A„ + l3{An+i + An-i) + 27r/3^„,i =0, n = 0, 1, 2, 



(3.1) 



(3.2) 



(3.3) 



(3.4a) 



(3.46) 



Proof. Equation (2.3) implies equation (3.2). Thus the problem reduces to solving an initial-boundary value 
problem for the heat equation (2.4)-(2.5). 
Let 



(p{x,n) = / ip{x,t)e~^"'*dt, 
Jo 



(3.5) 



^ OO if 

= ^ E 'Pix,n)e^-' = — U{x,0)+J2 

n=—oo \ 1 



(^(x,n)e*"*+v3(a;,-n)e-*"*) [. 



Equation (3.5) and reality imply 



ip{x, — n) = ip{x, n) 



Thus 

Equations (3.5) and (2.4a) imply 



= i- |^(x,0) + f;((^(a;,n)e-* +^(a;,n)e-'"*)| . 



'Pxx + 2c<^a; = 2 / 

Jo 

Thus using periodicity, we find 



27r 



(^t(a;,i)e-*"*rft = 2 \ip{x,t)e-'''^\^^ + irvf 



(pxx + 2c(^a; - 2m(^ = 0. 



(3.6) 



(3.7) 



n = 0: 

<^ = ci+^oe^"% Ao = -2c; 
the boundness requirement as a; ^ oo implies ci = 0. Thus 

n=l,2,---: 

.^(x, n) = .4„p/"^, A,^j + 2cA„ - 2in = 0. 
Substituting (3.8) into (3.6) we find (3.3). If the boimdary condition (2.2a) is vaUd, then 



(3.8a) 
(3.86) 



^„ = V>(0,n)= / ^{0,1)6 
Jo 



e-*"*dt, n = 0,1,2,- ••, 

Jo 

which is equation (3.4a). 

If the boundary condition (2.2b) is valid, then multiplyng this equation by e~*"* and integrating from to 27r we 
find 



e-"'' [2v5^(0, t) + /3(e" + e-")y^(0, t) + /3(e'* + e"^*)] =0, n = 0, 1, • • • 



2^.^(0, n) + I3[ip{0, n + 1) + <^(0, n - 1)] + 27r/3(5„,i = 0, n = 0, 1, • • • . 

<^x(0,n) = A„A„, 

the above equation becomes (3.4b). 



Thus 
Using 



4 The Boundary Value Problem (2.1)- (2.2a) with c = P'^k as P ^ 

Let 

c = l3'^k. (4.1) 

The solution of this problem is given by (3.2) and (3.3) where An is defined by (3.4a). The integral (3.4a) can be 
computed explicitly. For simplicity we only compute it as /3 ^ 0: 



/■27r 

Jo 



02 

Pcost+^{costf +0{l3^) 



\z\ = l 



2i 



1 / 1 1 J 

z -\ — ] az 



\z\ = l z 



n+l 



^(^+7) + ^(^' + 72+2)+0(/3^) 



dt, n = 0, 1, • • 

dz 



8i 



„^K + 2+i)rfz + 0(/3^). 



Thus 



Hence 



Also 



A„. = 2i7r 



TTT^n,! + (2(5„,o + (5n,2) 



+ 0(/3'). 



Ai=7r/3 + 0(/33), Ao = ^+0(/33), = !^ + 0(/33). 



Therefore equation (3.3) implies 



1 f/g^TT 



2tt 2 



+ 7r/3e" 



g-(l+i)x+it _j_ g-(l-i)x-it 



^ 4 



g-\/2(l+i)a;+2it _|_ g-%/2(l-i)a;-2it 



+ 0(/33 



6 



1. x = 0{l) 

r_(i + i)e-(i+i)a:+it _ (1 _ ^)g-(i-i)x-it] ^ |r-V2(l + i)e-v^(i+i)^+2it) _ ^(1 _ j)e-v^(i-i)x-2it] ^ o/^2) 

-71 ^ Z 

2 + /J[e-(i+i)^+'* + e-(i-»)^-'t] + 0(/32) 



or 



"=2 



(1 + *)e-(i+^)-+'* - (1 - i)e-(i-')--'* + ^ + ' + '''^ ' (4.2) 



_2e-2- -{l + i) e-2(i+i)-+2»t _ (1 _ i)e-2(i-i)-2it| + o{p^) 

2. x = 0{p-^) 

5 The Analysis of (2.1)-(2.2b) with c = P^k, k ^ 1/8, as ^ 

The solution of the hnear homogeneous equation (3.4b), with the requirement that 

An as n — > oo (5.1) 

is unique. This solution is given by 

An = -27TPFi + {-l)''AoFi---Fn, (5.2) 

where Fn is defined by 

^" = 2A„ / 1 — (5-3) 

3 

and the real constant Aq is determined from 

-4cAo + P{Ai + Ai) = 0. (5.4) 

Indeed, the definition F„ implies 



Solving this equation for 2A„//3 we find 



2A„ _ 1 



Substituting this expression into equation (3.4b) we obtain 

An+l+Fn+lAn + ^{An+FnAn-l) = -2'Kl35n,l, n = 0, 1, 2, • • • (5.5) 

Letting 

Gn = An + FnAn-U (5-6) 

equation (5.5) becomes 

Gn+i + = -27r/3(5„, 1. (5.7) 

Since Fn ~ 2a~' ^ cc, it follows that there does not exist a homogeneous solution of equation (5.7) that decays 
as n — > oo. Thus G„ = — 27r/3Fi(5„, 1, and equation (5.6) yields 

An + FnAn-1 = -27r/3Fi(5„,i. 
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The unique solution of this equation is given by (5.2). Evaluating equation (3.4b) at n = and using A_i = Ai we 
find (5.4). 

Equation (5.2) implies that 



n ^ 00, 



(5.5 



thus the series (3.3) converges. 



The asymptotic behavior as /3 — > 

The definition of A„ (equation (3.1)) together with the requirement that Re A„ < 0, imply that if c = /3^fc, then 



A„ = -^/^(l + ^) + 0(/32). 

Equations (5.2), (5.3) yield 

Ao = ao + 0(/3), Ai = /3ai + 0(/?'), ^2 = /3'a2 + 0(/3^), 
Equation (3.4b) with n = 0, 1, becomes 

n = : Ai+Ai- 4k(3Ao = 

n = 1 : 2AiAi + /3(A2 + Aq) + 27J-/3 = 0. 
Substituting equations (5.10) into (5.11), (5.12), we find 

oi + ai = Akao 

-2(1 + i)ai + ao + 27r = 0. 

Equation (5.13b) yields 



ai 



(ao + 27r)(l-i) 



Substituting this expression in equation (5.13a) and using that ao is real, we find 

ao 



8k -1' 

Substituting (5.10) into the expression for ip{x,t) (equation (3.3)) we obtain 



(5.9) 
(5.10) 

(5.11) 

(5.12) 

(5.13a) 
(5.136) 

(5.14a) 
(5.146) 



•^{x^t) = — 



1 r 27r 



27r 8fc - 1 



AkiT 
8k-l' 



-kfj^x 



1. x = 0(/3-2) 



2. x = 0{l) 



(1 - j)e-(i+')^+'* + (1 + i)e-(i-')^-**l p i + OiP^) 



2/3fc 



u{x,t) 



1 + (8A: - l)e2/3''=^ ■ 



g-(l+i)a:+it _j_ ^-{l-i)x-it 



+ 0{(3). 



(5.15) 



(5.16) 
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Appendix 

The Analysis of equation (3.4b) 

Let An satisfy the linear homogeneous difference equation 

2A„ 



Ar,+l + ^An + Ar,-l=0, n = 2,3, 



(Al) 



where /? is a constant and A„ — > oo monotonically as n ^ oo, for example A„ = A„, where A„ is defined by equation 
(3.1). We shah show that 

n—l r n—1 ^ n— 1 



An 



Cl 



n— 1 . 



C2 



n ^ 00. 



Indeed, we first make the change of variables 



Thus equation (A.l) becomes 



n-l 

An = Bnl[ 



T-> ^7") -^n—l 

-Dn+1 + -WT~^n + 







We write this equation in matrix form. 
The eigenvalues and the eigenvectors of the matrix Vn are 



1 2_ 



n = 2,3,- 








), 


) 




-t>n+l / 



Mr, 



± 



1 



0Kn V P^Al KK-i 



Using the gauge transformation 



we find 



Equation (A. 6a) implies 



1 1 



= T„$„, Tn = 
$„ = M„r-iT„_i$„_i, Mn 



IJ'n 



Mn = - 



Mn 



2AnA2_, 
2 



i + o(^)), 



l + 0(^))- 



These estimates and the explicit form of Tn imply that 



T-iT„_i=7 + 0(-l), 
thus the WKB approximation of equation (A. 8a) yields the two equations 

fn ~ Mn/"— 1' fn ~ 



{A.2) 



(A.3) 



(^.4) 



(^.5) 



(^.6) 



{A.7) 



(^.8) 



(A.9) 
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We recall that the general solution of the first order difference equation -F„+i = QuFn, is given by 

Fn= 9l - ■ ■ 9n-l- 

Thus to the leading order as n — > oo, 




Hence using (A.7) and (A.3) we find (A.2). 
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